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We consider a two-dimensional strongly localized system defined in a half-space and whose transfer 
integral in the edge can be different than in the bulk. We predict an unbinding transition, as the 
edge transfer integral is varied, from a phase where conduction paths are distributed across the 
bulk to a bound phase where propagation is mainly along the edge. At criticality the logarithm of 
the conductance follows the Pi Tracy-Widom distribution. We verify numerically these predictions 
for both the Anderson and the Nguyen, Spivak and Shklovskii models. We also check that for a 
half-space, i.e., when the edge transfer integral is equal to the bulk transfer integral, the distribution 
of the conductance is the Pi Tracy-Widom distribution. These findings are strong indications that 
random signs directed polymer models and their quantum extensions belong to the Kardar-Parisi- 
Zhang universality class. We have analyzed finite-size corrections at criticality and for a half-plane. 
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The probability distribution function of the conduc¬ 
tance g of quantum random hopping models, such as 
the Anderson model, has been much studied. In one 
dimension, it has been shown that all the cumulants of 
Ing scale linearly with system size [1]. Thus, the dis¬ 
tribution function of In g approaches a Gaussian form 
for asymptotically long systems and is fully character¬ 
ized by two parameters, the mean (In g) and the variance 
= (In^ g) — (In g)^. Both parameters are related to each 
other, supporting the extension of the single parameter 
scaling hypothesis [5] to the distribution function of the 
conductance [3]. 

In two-dimensional (2D) systems, calculations of the 
conductance distribution function are possible in the 
metallic regime, thanks to the non-linear sigma model 
[4], but not until now in the strongly localized regime. 
However, two of us have argued, and demonstrated nu¬ 
merically, that in that regime L ^ being the local¬ 
ization length, the conductance takes the form I3E] 
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where a is a constant and y a random variable with 
a Tracy-Widom (TW) cumulative distribution function 
(CDF). The TW distributions are CDF for the largest 
eigenvalues of large Gaussian random matrices [7] . It was 
found that the random variable y depends on the geome¬ 
try of the model. For narrow leads, the TW distribution 
associated to the Gaussian unitary ensemble y = X 2 is 
the CDF F 2 {x), i.e., Prob(y 2 < x) = F 2 {x) = F{xY, 
where F{x) is defined as F{x) = exp{—^ — 

x)u(s)^ds} and u(x) is the solution of the Painleve II 
equation u"{x) = 2u{x)^ + xu{x), with u{x) —>■ —Ai(x) 
as X —>■ -boo, Ai(x) being the Airy function. 

The chain of arguments leading to Q goes as follows. 
It was argued by Nguyen et al. (NSS) [8] that quan¬ 
tum interference effects in the localized regime are faith¬ 
fully described by retaining only the shortest or forward¬ 


scattering paths. Medina and Kardar [3] studied in detail 
the NSS model, focusing on its formulation as a model 
of directed polymers (DP) with non-positive Boltzmann 
weights. They found that the variance of the tunneling 
probability increases with distance as for 2D sys¬ 
tems. This suggests that in 2D the DP with non-positive 
Boltzmann weights should be in the same universality 
class as the standard (i.e. positive weight) DP problem, 
well known to exhibit a wandering exponent of 2/3. A 
qualitative replica argument led to the same conclusion 

la¬ 
in the mean-time, much progress happened in the 
study of growth models in the Kardar-Parisi-Zhang 
(KPZ) universality class [TU]. The standard DP prob¬ 
lem belongs to this class (the KPZ height h maps onto 
InZ where Z is the DP partition sum). Remarkably, it 
was found that the TW distribution arises at large time 
scales (i.e. long polymer length) in discrete models in the 
KPZ class. For example, the height in the polynuclear 
growth model m , the optimal energy in DP models 
and the length of the longest increasing subsequence in 
a random permutation m all follow TW distributions. 

This body of facts thus led to the conjecture 0 with 
the 1/3 exponent of the KPZ class, and to its numer¬ 
ical verihcation. More recently the KPZ equation and 
the DP problem have been solved directly in the contin¬ 
uum, using integrability by the Bethe Ansatz of an as¬ 
sociated quantum boson model [IMZl- Various bound¬ 
ary conditions where treated, suggesting new tests for 
the conjecture Q with various lead geometries. In par¬ 
ticular the continuum DP problem in a half-space was 
found [E] to lead to the F4 TW distribution, associated 
to the Gaussian simplectic ensemble, in agreement with 
earlier results for discrete models uninD]. This distri¬ 
bution verifies F^^x) = F{x){E{x)~^ E{x))/2 where 

A(x) = exp{i/+°' ’ u{s)ds}, F{x) and u{x) are the func¬ 
tions defined above. It was thus conjectured [18] that in 
a half-space the conductance near the edge should be of 
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the form Q with x = X 4 oi CDF given by F 4 . 

The aim of this paper is to study the logarithm of the 
conductance In g in 2D systems in the strongly localized 
regime in order to (i) verify numerically the conjecture 
that the CDF of for a half-plane is the F/i{x) function, 
and (ii) predict and analyze the transition from 2D to ID 
as the hopping amplitude t at the edge is increased, fa¬ 
voring propagation along the edge. Exactly at criticality 
we observe that the CDF of In g is the TW distribution 
^ 1 ( 3 ::) = F{x)E{x) of the Gaussian orthogonal ensem¬ 
ble. This transition is a generalization of the unbinding 
transition studied for positive weight DP [201 m]. Re¬ 
markably, the full crossover distribution obtained there 
fits our data for the whole parameter range, providing a 
very delicate test that the random sign DP belongs to 
the KPZ universality class. This is all the more precious 
since at present no integrable system has been found to 
solve the NSS model. 

We focus on the Anderson model on a square sample 
of finite size L x L described by the Hamiltonian 

H — ^ ^ T ^ ^ A h.c. , ( 2 j 

i {i,j) 

where the operator a\ (oi) creates (destroys) an electron 
at site i of an square lattice and is the energy of this 
site chosen at random G [—VF/2, bF/2], with W the 
strength of the disorder. The double sum runs over near¬ 
est neighbors. The hopping matrix element tij is taken 
equal to 1 everywhere, except between the sites along one 
edge where is equal t > 1 (see figure [^. The value in the 
bulk sets the energy scale, while the lattice constant sets 
the length scale. The unit of conductance is 2e^/h. 




FIG. 1: Scheme of the geometrical arrangements of the a) 
Anderson model, and the b) NSS model. Blue dots are site 
with random disorder energies and red dots sites on the leads, 
without disorder. Thin lines correspond to the bulk hopping 
strength equal 1, while thick lines along the edge to t. 

For the Anderson model we have calculated the con¬ 
ductance through Landauer’s formula in terms of the 
transmission between perfect leads, arranged as schemat¬ 
ically shown in Fig. This is obtained from the Green 
function, which can be calculated propagating layer by 
layer [22]. It is convenient to propagate along the di¬ 
rection perpendicular to the leads, starting from the op¬ 
posite edge, so that each calculation of the bulk Green 
function can be used for any value of t. 


The Green function between two sites a and b can be 
written in terms of the locator expansion 

HGi&) = En:^> (3) 

r ier * 

where the sum runs over all possible paths connecting 
the two sites a and b. The NSS model assumes that, in 
the strongly localized regime, is dominated by the 
forward-scattering paths and only takes into account the 
contributions of such paths to the transmission amplitude 
between opposite corners of a square lattice. As we want 
to study a half-space, we consider a triangular sample as 
represented in Fig. In order to mimic the Anderson 
model, we choose the site disorder energy at random in 
the interval G [—IF/2, IF/2], but if \E — ti\ < I, we set 
\/\E — ei\ = I, partially incorporating multiple scattering 
effects. We note that at least for the planar symmetry 
this choice of disorder does no change the universality of 
the problem [25] . 

The calculation of the quantum amplitude in the NSS 
model is formally similar to the calculation of the parti¬ 
tion function of DP in a random potential 

Z = E exp I -/3 E ^4 > (4) 

r I ier J 

where /3 = 1/kT, hi is a random site energy and F runs 
over all possible configurations of the DP. Equations ([^ 
and @ are equivalent provided that we can identify —jShi 
with ln(£’ — Ci). While in DP the disorder energies hi are 
real, in the quantum case E — Ci does not have to be 
positive, implying a more general problem. 

It is interesting to study first the conductance distri¬ 
bution of the NSS and Anderson models in a half-space 
and check whether the conjecture of Ref. [15], based on 
the exact result for the continuum DP, that the GDF of 
X is the TW function F 4 is verified by our two models. 
In Pig. [^ we plot histograms of h\g for the NSS model 
as a function of z = {hig — A)/B, where A and B are 
chosen so that z has the same mean and variance as E 4 . 
The lateral dimensions are L = 2500 (blue dots), 5000 
(red dots) and 10000 (black dots), and the number of re¬ 
alizations is 4 X 10®. The solid line corresponds to the 
TW function E'^. We see a perfect agreement between 
our numerical results and the TW distribution for more 
than four orders of magnitude. A similar agreement is 
found for the Anderson model (not shown), taking into 
account that it corresponds to smaller system sizes. 

The inset of Fig. [^ we plot the skewness of In g ver¬ 
sus which according to the discussion below is 

expected to be the leading order correction. The red line 
corresponds to the NSS model, while the dots to the An¬ 
derson model (IF = 10 blue, 20 magenta). The dashed 
line is Sk 4 = 0.165509... The skewness of the numerical 
results for both models tends to the predicted value. 
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FIG. 2: Histogram of Ing versus the scaled variable 2 for 
three sizes and a disorder VF = 10 of the NSS model. The 
continuous line corresponds to F^{z). Inset: skewness vs. 

for the NSS model (red curve) for the Anderson model 
with IF = 10 (blue dots) and 20 (magenta dots); the dashed 
line is Sk 4 . 


We now argue theoretically that the NSS and the An¬ 
derson models exhibit a phase transition between (i) a 
2D phase where conduction paths wander unboundedly 
in the half-space and the fluctuations of In 5 at large L 
are described by 0 with X = X 4 i ^nd (ii) a quasi-ID lo¬ 
calized phase where the conduction paths have ID char¬ 
acter at large scale, and In 5 has a log-normal distribu¬ 
tion with fluctuations scaling as i.e. 0 still holds 

but with a size dependent variable X ^ XL, scaling now 
as XL L~^l^ . Exactly at the transition t = the 
fluctuations of In g at large L are described by 0 with 
now X = Xi distributed according to Fi. Hence the 1/3 
KPZ exponent holds up to and at the transition. This 
transition can be described as an unbinding transition 
of conduction paths. In the standard DP the unbinding 
transition was predicted in and worked out in detail 
in the framework of symmetrized random permutations 
in [50] • Before using these results below, let us first give 
the physical picture and the scaling arguments. 

In the unbound phase, t > tc, the paths wander in 
the 2D space, but come back from time to time to 
the boundary (they still feel the boundary hence the 
change from Fi to F 4 ). Their typical transverse wan¬ 
dering is Ins the bound phase, the paths return 

to the wall, with a typical length which diverges at 
the transition as ^ {t — Near the transition 

there are thus independent pieces of paths, each 

fluctuating as hence the variance should behave as 
(In 5 ) 2 '' = ~ {t-tc)L. Far from the 

transition this picture predicts that the skewness should 
decay as Sk~ (L/^„)“^/^, since in the ID phase the cen¬ 
tral limit theorem holds for Ing and all cumulants scale 
as L i.e. (Ing)” ^ L. 

Let us now consider the critical scaling around the 


transition point |20] . Our prediction is that Ing always 
behaves as in Q with the random variable x having 
a CDF Prob(x < x) = F{x,w = —cw) that depends 
on the scaling variable w = {t — tc)L^^^, c > 0 being 
an (unknown) proportionality constant. The universal 
crossover function F{x, w) (corresponding to F^(x, w/2) 
in Ref. 1201) satisfies 

F{x, w) = ]^F{x) {{f{x, w) - g{x, w))E{x)~'^ (5) 

F{f{x,w)+g{x,w))E{x)) 

where the functions f{x,w) and g{x,w) are the ones of 
[51] and verify: 

d,^f{x,w) = u{x)'^f(x,w)-{wu{x)+u'{x))g{x,w) 

dwg{x,w) = —{wu{x) — u'{x))f{x,w) + 

{w'^ — X — u^{x))g{x,w) (6) 

and are subjected to the initial condition at criticality 
f{x, 0) = —g(x, 0) = E(x)'^. Then, the CDF at criticality 
is J^(a;,0) = Fi{x). In the limit w —>■ -boo, one has 
/ —>■ 1 and g —t 0, hence F{x, 00 ) = F 4 . In the opposite 
limit, w —>■ —00 [51] F{w‘^ + y,y\w\,w) —)■ erf(g) and one 
recovers Gaussian fluctuations for Ing. 

To study the transition and to characterize the two 
phases, we analyze the skewness of the distribution of In g 
for the NSS model as we vary t. The skewness should 
tend to 0 in the localized phase, to Ski = 0.293464... 
at the transition and to Sk 4 in the unbound phase. In 
Fig. [3] we plot the skewness as a function of size on a 
logarithmic scale for several values of the edge strength 
t. The upper horizontal line corresponds to the expected 
value of the skewness at the critical point. Ski, and the 
lower horizontal line to the limiting value in the bulk 
phase, Sk 4 . The black curve represents the skewness at 
the critical value G ~ 1.613(1) and tends to the upper 
horizontal line. Ski. Solid curves are for t > F and after 
reaching a maximum start decreasing, eventually tending 
to zero. Dashed lines correspond to t < tc and tend to 
Sk 4 , middle horizontal line in Fig.[^ 

It is possible to scale the raw data for the skewness at 
large values of L into a universal curve using as scaling 
variable w = (t — tc)L^^^. Looking at the behavior of 
the skewness for the critical value tc in Fig. [^ it is clear 
that finite size effects are important. To take them into 
account we assume Sk(t, L) « Sk{w)f{L), where f{L) 
incorporates finite size effects in a simple way f{L) = 
Sk(tc, L)/Ski. In Fig. |^we plot this renormalized values 
Sk(?Zi) for several values of t. All curves are plotted for 
the range 200 < L < 10"^. The black dot represents the 
transition point and so is equal Ski and is placed at the 
vertical axis, t = tc- The dashed line is the theoretical 
prediction for the scaling function, i.e., the solution of 
Eqs. (|^, using a value c = 0.9 for the multiplicative 
constant. It presents a maximum in the localized phase 
well reproduced by our simulations. The blue dotted 
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FIG. 3: Skewness versus L on a logarithmic scale for the val¬ 
ues of t shown in the figure. The red dashed line corresponds 
to Ski and the dashed black horizontal line to Ski. 

line on the right represents the limiting behavior in the 
bound phase, with the fitting parameter b = 

0.55. The overlap and the agreement with theoretical 
expectations is quite good, specially noting that the only 
free parameter to overlap all the curves is tc- 


In g (here) and of the interface height h (there, time being 
denoted as L) and the scaled cumulants /c„ = 
which converge to constants. The finite size corrections 
of the kn where found (there in the bulk) to scale as 
for n = 1 and as at most for all n > 2 [^[5^. We 

study these corrections for the NSS model. In the top 
panels of Fig. we plot ^2 (left) and k^ (right) versus 
2 ^-2/3 Jqj, half-space, t = 1. We note that ks follows 
a very good linear behavior, while ^2 shows some curva¬ 
ture. Based on the curvature of /c 2 it is difhcult to de¬ 
termine its leading finite size corrections exponent, but 
the exponent —2/3 produces an extrapolated value for 
the skewness (Sk= 0.166) very close the the theoretical 
expectation, indicating that finite size corrections for the 
half-space are similar to those observed for bulk growth 
models in KPZ class. In the lower panels of Fig.[^we plot 
^2 (left) and k^ (right) at the critical point versus 
The smooth behavior of both curves indicate that finite 
size corrections exhibit a distinct and novel behavior at 
criticality. From the intercepts with the vertical axis of 
linear fits at large L, we obtain Sk= 0.297, confirming 
that the CDF is the TW function J^i. The histograms 
of In 5 at the critical point agree fairly well with these 
theoretical predictions. 



FIG. 4: Skewness scaled as a function oi w = (t — tc)L^^^ for 
several values of t. The black dot corresponds to Ski at the 
transition and the horizontal dashed line to Ski. The dashed 
curve is the theoretical prediction for the scaling function and 
the dotted curve the predicted asymptotic behavior at large 
w. Inset: skewness as a function of t for the Anderson model. 

The results for the Anderson model confirm all predic¬ 
tions, although the maximum size that can be calculated, 
L = 400, is still relatively far from the asymptotic behav¬ 
ior. In the inset of Fig. we plot the skewness as a func¬ 
tion of t for several system sizes of the Anderson model 
and a disorder IT = 20. A peak around t « 1.8 develops 
with size, but its maximum is still far from Ski. This 
results are fully consistent with NSS for similar sizes. 

It is interesting to compare the finite size corrections in 
our models with those for discrete growth models and ex¬ 
periments in the KPZ class. Consider Kn the cumulant of 


FIG. 5: Leading order corrections of the second fc 2 (left) and 
third fcs (right) cumulants in the NSS model for a half-space 
versus (top panels) and at the unbinding transition 

versus (bottom panels). 


We have shown that the NSS and the Anderson mod¬ 
els in a half-plane undergo an unbinding transition as the 
hopping amplitude at the edge is varied. We show evi¬ 
dence that the analytical expressions for the conductance 
distribution functions and the scaling functions at the 
unbinding transitions in localized two-dimensional quan¬ 
tum systems are obtained from Eqs. (5][6). We have also 
studied the corresponding finite-size corrections. Similar 
unbinding transitions can occur in 3D systems when con¬ 
duction through a surface or a line is favored. They may 
be relevant to the behavior of edge states. 
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